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ABSTRACT 

We explore the possibility that neutron stars accreting from the winds of main-sequence stellar com- 
panions account for a significant fraction of low-luminosity, hard X-ray sources (Lx < 10 35 ergss _1 ; 
1-10 keV) in the Galaxy. This work was motivated by recent Chandra observations of the Galactic cen- 
ter by Wang et al. (2002). Our calculations indicate that many of the discrete X-ray sources detected in 
this survey may be wind-accreting neutron stars, and that many more may be discovered with deeper 
X-ray observations. We propose that an infrared observing campaign be undertaken to search for the 
stellar counterparts of these X-ray sources. 

Subject headings: stars: neutron — X-rays: stars 
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1. INTRODUCTION 

Of all the binary X-ray sources in the Galaxy that con- 
tain a neutron star (NS), the most abundant, though typ- 
ically not the most luminous, are systems in which the NS 
accretes from the wind of a main-sequence stellar compan- 
ion. Immediately following the supernova (SN) explosion 
that accompanies the formation of the NS, the companion 
is, in most cases, relatively unevolved. Stellar masses rang- 
ing from a few to several tens of solar masses are likely, 
with corresponding lifetimes of ~ 10 7 -10 8 yr. The long 
phase of wind accretion onto the NS is a fairly quiet pre- 
lude to the more dramatic Roche-lobe overflow that fol- 
lows, where the NS may appear as a low/intermediate- 
mass X-ray binary (L/IMXB), or, more commonly, as a 
short-lived high-mass X-ray binary (HMXB) before being 
engulfed by its companion. 

Well known examples of wind-accreting NSs (WNSs) 
are present in HMXBs, which are typically assumed to 
contain > 10- Mq stars. It has long been suggested (e.g., 
Rappaport & van den Heuvel 1982) that there is a large 
unobserved Galactic population of HMXBs with steady X- 
ray luminosities of Lx ^ 10 35 ergss _1 . Most importantly 
for the present work, WNSs in HMXBs have been pro- 
posed (e.g., Ogelman & Swank 1974; Koyama et al. 1989) 
to account for much of the so-called "Galactic ridge" of 
previously unresolved, hard (~ 1-10 keV) X-ray emission. 

The nature and contributions of diffuse and discrete X- 
ray sources to the X-ray spectrum and total X-ray lumi- 
nosity of the Galactic ridge has been a point of debate 
for decades. Prior to Chandra, X-ray satellites that have 
surveyed the Galactic plane (e.g., ASCA, RXTE) lacked 
the angular resolution and/or sensitivity to identify point 
sources near the Galactic center. Wang, Gotthelf, & Lang 
(2002; hereafter, WGL02) recently reported on a high- 
resolution survey of a 0.8° x 2° field about Sgr A*, carried 
out with Chandra/ ACIST in the energy range of 1-8 keV. 
Throughout most of the observed field, the sensitivity for 
detecting discrete sources is <Sx,min ~ 10~ 13 ergs s _1 cm~ 2 , 
corresponding to Lx ~ 10 33 ergss _1 at the distance (~ 
8.5 kpc) of the Galactic center. 



WGL02 estimate > 500 previously undetected point 
sources in the Chandra survey, the nature of which are un- 
known. Sources with X-ray emission in the 1-3 keV band 
are probably within several kiloparsecs of the Sun. Many 
of the sources with energies > 3keV, for which the softer 
X-rays have been absorbed by the interstellar medium, are 
nearer to or beyond the Galactic center. WGL02 suggest 
that as many as ~ 50% of the hard point sources may 
be background active galactic nuclei (AGN), although the 
true fraction may be closer to ~ 10-20%, as estimated 
from the Chandra Deep Field data (Brandt et al. 2001). 
This leaves > 100 Galactic objects that are most likely a 
population of binary X-ray sources, including black- hole 
binaries (BHBs), LMXBs, cataclysmic variables (CVs), 
and WNSs. In the next section, we devote a short dis- 
cussion to each candidate class of X-ray binary, and we 
support the case that most of the hard (> 3keV) Galactic 
point sources detected by Chandra are WNSs. 

As we argue below, the field surveyed by Chandra en- 
compasses roughly 1% of the stars in the Galactic disk. 
This nontrivial fraction may translate to > 100 WNSs 
(see § 3), and thus the survey may provide an important 
partial census of these objects. However, even if future ob- 
servations reveal that few of the hard sources detected by 
Chandra are WNSs, the main ideas of this paper provide a 
foundation for constraining their numbers and properties. 

2. THE NATURE OF THE POINT SOURCES 

The total numbers for each class of X-ray binary in 
the solid angle of 1.6 deg 2 observed by Chandra may be 
estimated as follows. We suppose that the space den- 
sity of stars in the Galactic disk is given by n(R, z) oc 
exp(— R/Rq) exp(— \z\/zo), where R is the galactocentric 
radius, and z is the displacement from the Galactic mid- 
plane. The scale radius, Ro, is taken to be 4 kpc (van der 
Kruit 1987), and the vertical scale height is zq > 100 pc. 
Upon integrating n{R, z) over 1.6 deg through the Galac- 
tic center, we find that the field contains < 1% of the 
Galactic population of hypothetical X-ray sources. Thus, 
the detection by Chandra of > 100 Galactic sources im- 
plies a total number of > 10 4 in the entire Galaxy. We 
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now consider each of the candidate X-ray binaries in turn. 

Various semi-empirical and theoretical estimates give to- 
tal numbers of LMXBs and BHBs (quiescent and active) in 
the Galaxy at ~ 10 3 each (e.g., Verbunt & van den Heuvel 
1995; Romani 1998). Therefore, it seems that at most of 
order 10 LMXBs and BHBs could have been detected in 
the WGL02 survey. This estimate is consistent with the 
< 20 LMXBs and one BHB observed in the surveyed field 
(see Liu, van Paradijs, & van den Heuvel 2001). 

There are probably ~ 10 6 CVs in the Galaxy (e.g., How- 
ell, Nelson, & Rappaport 2001), and therefore ~ 10 4 CVs 
in the field surveyed by WGL02. Most observed CVs have 
Lx ^ 10 32 ergss _1 in the bandpass of 0.1-2.5keV, and 
decreasing power toward higher energies. Luminous CVs, 
with Lx ~ 10 33 -10 34 ergs s -1 , which should comprise only 
~ 1% of the intrinsic population (e.g., Howell, Nelson, & 
Rappaport 2001), likewise have relatively soft X-ray spec- 
tra. Therefore, we expect that few CVs are detectable be- 
yond several kiloparsecs for Sx.min = 10~ 13 ergs s _1 cm~ 2 
(1-8 keV), because of their low total Lx, as well as the 
heavy interstellar absorption for energies < 3 keV. Inte- 
grating n(R, z) over the 1.6 deg surveyed by Chandra, to a 
distance of D — 3 kpc from the Sun, we calculate that sev- 
eral tens of CVs, and perhaps of order 100, may have been 
detected. Therefore, CVs may make up a sizable fraction 
of the 1-3-keV sources in the WGL02 mosaic image. 

Arguably, the most compelling hypothesis for the hard 
(> 3keV) Galactic X-ray sources detected by Chandra 
is that they are WNSs. Observed HMXBs exhibit a 
wide range of X-ray luminosities (~ 10 33 ~ 10 38 ergss _1 ) 
and hard, nonthermal X-ray spectra that extend beyond 
lOkeV (Nagase 1989). Many (~ 40%) of these systems are 
transient, with short (less than several months) outburst 
phases where Lx > 10 36 ergss _1 , and recurrence times 
that are often years or decades. During its essentially qui- 
escent state, a transient HMXB may have Lx < 10 35 , due 
to steady accretion from the stellar wind by the NS, such 
as in the case of the Be/X-ray binary X Per/4U 0352+30 
(Dclgado-Martf et al. 2001). Theoretical calculations (e.g., 
Meurs & van den Heuvel 1989) indicate that tens of thou- 
sands of such low-luminosity WNSs may currently inhabit 
the Galaxy, implying that hundreds of WNSs may popu- 
late the 0.8° x 2° field about the Galactic center. In the 
next section, we use binary population synthesis to esti- 
mate the number of potentially observable WNSs. 

3. THE X-RAY FLUX DISTRIBUTION 

An estimate of the number of WNSs with Sx > <Sx,min 
in the field surveyed by Chandra requires that we calcu- 
late their X-ray flux distribution. The X-ray luminosity 
of an individual source depends upon the parameters that 
characterize both the stellar wind and the binary system 
that contains the WNS. We have conducted a Monte Carlo 
population synthesis study of the formation of binaries 
that consist of a NS and a stellar companion. The three 
main steps of our population synthesis calculation are enu- 
merated below. More detailed descriptions of massive bi- 
nary stellar evolution and the elements of our population 
synthesis code are given in Pfahl, Rappaport, & Podsiad- 
lowski (2001; hereafter, PRP01) and Pfahl et al. (2001). 

1) In the Galactic disk, each binary containing a NS de- 
scends from a massive primordial binary, where the ini- 



tially more massive and less massive stars are hereafter re- 
ferred to as the primary and secondary, respectively. We 
take the Galactic formation rate of massive binaries to be 
comparable to the core-collapse SN rate, TZsn ~ 10~ 2 yr -1 
(Cappellaro, Evans, & Turatto 1999). The initial pri- 
mary mass, taken to be Mu > 8M Q , is chosen from a 
power-law IMF, p(M u ) oc M^ 2 5 . The initial secondary 
mass, M-2i, is chosen from a flat distribution of mass ra- 
tios, qi = M^i/Mu < 1. For simplicity, we assume that 
the primordial binaries have circular orbits (see PRP01), 
and we choose the initial orbital separation, a i: from a 
distribution that is flat in log a, . 

2) If the orbit is sufficiently compact (a, < 5 — 10 AU) 
that the primary evolves to fill its Roche lobe, we use 
analytic formulae (PRP01) to compute the orbital sep- 
aration following the subsequent phase of mass transfer. 
The mass ratio and the evolutionary state of the primary 
at the onset of Roche-lobe overflow are used to deter- 
mine whether the mass transfer is stable or dynamically 
unstable. Given some critical mass ratio, which we take 
to be q c = 0.5, the mass transfer is assumed to be sta- 
ble if qi > q c and the envelope of the primary is mostly 
radiative when mass transfer begins, and dynamically un- 
stable if qi < q c or the primary has a convective enve- 
lope. We assume that the entire hydrogen-rich envelope 
of the primary is removed during mass transfer, whether 
stable or dynamically unstable, leaving only the primary's 
hydrogen-depleted core. For stable mass transfer, we sup- 
pose that the secondary accretes a fraction, (3 = 0.75, of 
material donated by the primary, and that the remain- 
ing mass escapes the system with a specific angular mo- 
mentum that is a = 1.5 times the orbital angular mo- 
mentum per unit reduced mass. The orbital separation 
increases or decreases during stable mass transfer by a 
modest factor of < 5-10 for reasonable values of a and (3 
(see PRP01). We assume in all cases that the secondary 
is "rejuvenated" due to the accretion, so that it emerges 
following mass transfer on the ZAMS appropriate for its 
new mass. The minimum secondary mass following sta- 
ble mass transfer is ~ q c 8 M Q + (36 M Q = 8.5 M Q , where 
we have used our reference values of q c and (3, and 8 M 
and 6 Mq are, respectively, our chosen minimum value 
of M\i and the corresponding envelope mass. Dynami- 
cally unstable mass transfer is accompanied by a common- 
envelope phase, wherein the secondary experiences a drag 
that causes the binary orbit to shrink in a time of < 10 3 yr. 
A fraction, r/QE ^ 1, of the initial orbital energy is avail- 
able to unbind the common envelope from the system. If 
insufficient energy is available, the two stars will merge. 
Otherwise, the envelope of the primary is dispersed, and 
the secondary emerges near the ZAMS, without having 
accreted any mass. The orbital separation may be ~ 100 
times smaller for binaries that survive dynamically unsta- 
ble mass transfer. A merger results in nearly every case 
where qi < q c and the primary's envelope is radiative when 
mass transfer starts. This implies that almost all systems 
that survive dynamically unstable mass transfer must have 
had initial orbital separations wide enough for the primary 
to grow to become a convective red supergiant. 

3) Upon exhausting its remaining nuclear fuel, the exposed 
core of the primary explodes as a Type lb or Ic SN. The 
impulsive mass loss and possibly large "kick" to the NS 
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strongly perturb, and may unbind, the binary. We assume 
that both the mass loss and the kick are instantaneous, and 
that the orientations of the kicks are distributed isotropi- 
cally. Two scenarios are considered for the distribution of 
kick speeds. Observations of isolated radio pulsars indicate 
that the mean NS kick speed may be > 100-300 km s _1 
(e.g., Hansen & Phinney 1997; Arzoumanian, Chernoff, & 
Cordes 2001). In our first kick scenario (hereafter, Kl), 
we apply a Maxwellian distribution in kick speeds with 
a mean of ~ 300 km s -1 to all NSs. Our second scenario 
(K2) was developed (Pfahl et al. 2001) to account for a new 
class of HMXBs with long orbital periods (P or b > 30 d) 
and low eccentricities (e < 0.2). For NS progenitors that 
are able to evolve into red supergiants (i.e., single stars 
and those in wide binaries), the kick speeds are drawn 
from a Maxwellian with a mean of ~ 300kms _1 , as in 
Kl. Therefore, our arguments in the last paragraph imply 
that for post-SN binaries where M 2 < 8M Q , the NS has 
received the conventional large kick. If the NS progenitor 
is in a binary system, and its envelope is removed while it 
is mostly radiative, which includes all systems where the 
mass transfer is stable, we utilize a Maxwellian distribu- 
tion with a much lower mean of ~ 30 km s . Kick scenario 
K2 yields a much larger number of post-SN binaries con- 
taining massive companion stars (> 8M©), as compared 
to the number of surviving systems with M 2 < 8 M . 

The relevant output of the population synthesis code 
is the post-SN secondary mass, M 2 , semimajor axis, a, 
and eccentricity, e, for each binary. We neglect tidal cir- 
cularization during the main-sequence evolution of the 
secondary, although this is important for systems with 
a (l — e ) ~ 30-40 Rq. Our results will not be greatly mod- 
ified if we include tidal circularization. 

Each WNS accretes from the wind of a relatively un- 
evolved stellar companion. Winds from early-type stars 
with masses of ~ 3-20 M are characterized by high 
speeds, v w ~ 1000 km s -1 , and low to moderate mass- 
loss rates, M w ~ 10" n -10~ 7 M yr" 1 . Kudritzki & Puis 
(2000) give the asymptotic wind speeds for stars near the 
main sequence, which vary from ~ 1.5 times the sur- 
face escape speed, v e , for stellar masses of ~ 3-10 M , 
and ~ 2.5v e for hotter, more massive stars. The wind 
speeds are quite uncertain, however, and we thus con- 
sider two simple cases in our simulations, namely v w = v e 
and v w — 2v ei independent of M 2 . For M w , we utilize 
the fitting formula of Nieuwenhuijzen & de Jager (1990): 
M w cx M 2 116 Lp 4 i?^- 81 , where L 2 and R 2 are, respectively, 
the luminosity and radius of the secondary. For both M w 
and v e , we substitute the time-averaged values of L 2 and 
i? 2 during the main-sequence phase — roughly twice their 
ZAMS values for a wide range in M 2 . 

We apply the standard Bondi-Hoyle-Lyttleton accretion 
scenario (Hoyle & Lyttleton 1941; Bondi & Hoyle 1944) to 
obtain the X-ray luminosity of each WNS, specialized to 
the case where v w is much larger than the orbital speed 
of the NS. We further assume that the wind is steady 
and spherically symmetric, so that the density varies as 
M^r -2 ^" 1 , where r is the instantaneous orbital separa- 
tion. The orbital time-averaged X-ray luminosity then 
scales as (L x ) oc eM m (i _2 u„ 4 (l - e 2 ) -1 / 2 . Here, e < 1 
is the efficiency for converting gravitational energy into 



hard (1-10 keV) X-radiation. For e = 0, a = 0.5 AU, 
v w = 1000 kms" 1 , and M w = lO~ 8 M yr- 1 , we find, af- 
ter including the multiplicative factors appropriate for an 
accreting NS, that (£x) ~ e • 10 33 ergss _1 . In the present 
work, we do not take into account the spin history of the 
NS and the centrifugal inhibition of accretion (e.g., Stella, 
White, & Rosner 1986), although this should be incorpo- 
rated into a more detailed study. 

The X-ray flux distribution of the WNSs is obtained 
by convolving the distance and X-ray luminosity distri- 
butions, which we assume are independent in this study. 
Since the probability of observing a certain WNS should 
be roughly proportional to the main-sequence lifetime, 
tms(-^2) ; of the secondary, we compute the distribution 
of intrinsic X-ray luminosities by accumulating the values 
of tms(-^2) for each bin in (Lx)- Using the same expo- 
nential form given in § 2 for the space density, n(R,z), 
of WNSs, we obtain the distribution of distances, D, from 
the Sun, for sources located within a solid angle of 1.6 deg 2 
centered on the Galactic center. We adopt fixed values of 
i?o = 4 kpc and zo = 200 pc as being typical for the post- 
SN binaries. In reality, the effects of SN mass loss and NS 
kicks are such that, for a given kick speed, binaries with 
lower-mass secondaries will receive larger systemic veloc- 
ities and reach greater distances from the Galactic mid- 
plane. A more detailed investigation should include the 
dynamical evolution of post-SN binaries in the Galaxy. 

Our binary population synthesis calculation yields the 
formation efficiency, Ffe, for post-SN binaries with M 2 > 
3Mq (a somewhat arbitrary minimum mass), where the 
remainder of the massive binaries have ultimately merged 
or have been disrupted due to the SN. The secondary stars 
in the surviving systems have a mean main-sequence life- 
time of (tms) ; so that the total number of WNSs in the 
Galaxy is ~ (tms^SN-Tte- Multiplying this total Galac- 
tic number by the fraction, J-"q, of systems encompassed 
by the observed solid angle, we obtain the intrinsic num- 
ber of WNS in the surveyed field. If a fraction, JT X , of the 
WNSs have Sx > Sx.min, as determined from the X-ray 
flux distribution, Chandra should have detected roughly 
{tms^Sn^fe^q^x such sources. 

Here, we define intermediate- and high-mass post-SN 
binaries by M 2 = 3-8 M and M 2 > 8M , respectively. 
When we apply the conventional kick scenario, Kl, we ob- 
tain J-pE ~ 0.05. The general importance of intermediate- 
mass secondaries is indicated by the rather large mean 
lifetime of (tms) ~ 50Myr for the post-SN binaries. The 
formation efficiency in kick scenario K2 is .T-fe ~ 0.2, and 
the mean lifetime is reduced to (tms) ~ 30Myr. We find 
that Tii ~ 0.01, so that the total number of WNSs in the 
solid angle surveyed by Chandra is ~ 250 and ~ 600 for 
kick scenarios Kl and K2, respectively. 

Examples of the X-ray flux distribution for the WNSs 
are shown in Fig. 1. The two curves in the top panel 
correspond to our two choices for the wind speed, both 
for kick scenario Kl. Kick scenario K2 does not yield an 
appreciably different flux distribution, aside from the nor- 
malization factor J-fe- The middle and bottom panels of 
Fig. 1 show the contributions of intermediate- and high- 
mass systems, which change considerably in the two kick 
scenarios for reasons given earlier. The fraction of WNSs 
with Sx > fix, min depends very much on the prescription 
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Fig. 1. — The X-ray flux distributions for a set of model as- 
sumptions. The top panel shows the distributions for v w = v e 
(right curve) and v w = 2v e (left curve), where kick scenario Kl was 
used. The middle panel shows the contributions of intermediate- 
mass (short- dashed) and high-mass (dotted) systems, for v e = 2v w 
and kick scenario Kl. The bottom panel is similar to the middle 
panel, but for kick scenario K2. 

for the stellar wind speed. For v w — v e and v w = 2v e , 
we find that .Fx ~ 0.4 and ~ 0.05, respectively, giving de- 
tected numbers of WNSs between ten and several hundred. 

Figure 1 shows that WNSs with intermediate-mass com- 
panions may make a dominant contribution to the flux dis- 
tribution over a wide range in Sx- This is largely due to 
the longer lifetimes and and lower wind speeds for lower 
secondary masses. While only a few of the observed so- 
called HMXBs (Liu, van Paradijs, & van den Heuvel 2000) 
have the mid- to late-B spectral types consistent with un- 
evolved intermediate-mass stars, it is certainly plausible 
that a much larger number of the known systems harbor 
such companions; in most cases the spectral subtype is not 
well constrained. Interestingly, the small fraction of these 
systems with M 2 < 4 M Q and P or b ~ 1-10 d will appear as 
iiviXBs (Podsiadlowski, Rappaport, & Pfahl 2002) when 



the secondary fills its Roche lobe and stably transfers mat- 
ter to the NS. For longer periods and larger masses, the 
mass transfer is dynamically unstable. However, as es- 
sentially all of the NS binaries with < 8 M Q form in 
much the same way, the entire population of WNSs with 
intermediate-mass companions may provide unique statis- 
tical information regarding the formation of IMXBs. 

4. DISCUSSION 

Chandra has detected ~ 1000 X-ray sources about the 
Galactic center. The next step is to determine the na- 
ture of these sources. Based on the discussions and results 
presented in the last two sections, we are in a position 
to make suggestions regarding future observational work, 
especially as it pertains to WNSs. 

A significant fraction — perhaps the majority — of the 
3-8-keV point sources detected by Chandra may be WNSs 
located within several kiloparsecs of the Galactic center. 
We propose that an infrared observing campaign be un- 
dertaken to search for stellar counterparts. Our models 
indicate that the majority of WNSs have companions with 
mid-0 to late-B spectral types, with unreddened if-band 
magnitudes of ~ 11-16 at the distance of the Galactic 
center, assuming these stars are near the ZAMS. The in- 
trinsic JHK colors of hot stars are nearly degenerate, and 
determination of their spectral subtypes thus requires a 
spectral-line classification system (e.g., Hanson, Conti, & 
Reike 1996). Purely photometric observations are an im- 
portant first step, however, as they should allow one to 
distinguish between stars and background AGN, even with 
the effects of interstellar extinction. 

We also suggest that it would be worthwhile to extend 
the Chandra survey in both sensitivity and angular cover- 
age. The X-ray flux distributions shown in Fig. 1 indicate 
that a factor of 10 increase in sensitivity may gain a factor 
of ~ 5-10 in the number of detected WNSs; of course, the 
number of AGN may likewise increase by a factor of sev- 
eral tens. Furthermore, extending the survey to Galactic 
latitudes |6| > 2° may reveal a gradient in the density of 
point sources, and thus give the fraction that are Galactic. 

We are grateful to F. Baganoff, N. Brandt, D. 
Chakrabarty, A. Juett, M. Muno, and Q. D. Wang for 
useful discussions. This work was supported in part by 
NASA ATP grant NAG5-8368. 
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